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Abstract. Recent advances in Fluid Dynamical modeling of heavy ion collisions are 
presented, with particular attention to mesoscopic systems, QGP formation in the pre 
FD regime and QGP hadronization coinciding with the final freeze-out. 
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1. Fluid Dynamics, and Local Equilibrium 

Fluid dynamics (FD) is probably the most frequently used model to describe heavy ion 
collisions. It assumes local equilibrium, i.e. the existence of an Equation of State (EoS), 
relatively short range interactions and conservation of energy and momentum as well as of 
conserved charge(s). Thus, it is a widely usable model. 

It can be derived most simply from the Boltzmann transport equation where, the con- 
servation laws, assuming local (approximate) kinetic equilibrium yield the equations of 
(viscous) fluid dynamics. If the local momentum distribution deviates strongly from local 
kinetic equilibrium the fluid dynamical approach (at least the one fluid one) is not appli- 
cable. This makes many believe that fluid dynamics has less applicability than transport 
models, like molecular dynamics models. 

One tends to forget other assumptions in transport models, i.e. dilute systems with 
binary collisions only, and consequently binary collisions. These, constraints limit the 
applicability of transport models, for example phase transitions (which include strong cor- 
relations, dense systems, and not only binary collisions) can hardly be described correctly 
with transport models. 



xxxxx-xxx/2000 
©2000 EP Systema Bt., Debrecen 



2 



L. Csernai et al. 



Classical FD models incorporate phase transitions in a trivial way, as the EoS, is given 
for a phase transition. Even more so, the FD approach can describe systems out of phase 
equilibrium, supplemented with a dynamical equation describing the dynamics of the phase 
transition, as local kinetic equilibrium in each phase is sufficient to apply FD.[[l|, ^ 

Phase transition dynamics is an involved question, even in macroscopic systems. First 
of all, phase transitions can be different. This expression may include, slow burning or 
deflagration, detonation, condensation, evaporation, and many other forms of transition. 
The basic conditions of all these transitions have, nevertheless, some similarities. These 
arise from the basic conservation laws and from the requirement of equilibrium. 

If we consider macroscopic stationary systems, asymptotically both the initial and final 
states are in local mechanical (pressure), thermal and chemical (or phase) equilibrium. The 
spatial extent and time-span of the transitional region depends, on the other hand, from 
many features of the transition. 

Most explicit dynamical calculations are performed for the homogeneous nucleation 
geometry as this is usually the mechanism which starts the transition and which is the 
slowest of all. 

In a dynamical situation the approach using the EoS including a first order phase tran- 
sition is identical both in situations involving compression or expansion. If the compression 
is supersonic, Shockwaves or detonation waves are formed, where the final new phase is 
immediately formed. The phase transition speed influences only the width of the shock 
front, but for slow dynamics and rapid phase transition the shock front width is primarily 
determined by the transport coefficients, viscosity and heat conductivity, and not by the 
phase transition speed. 

If phase transitions occur in small finite systems other dynamical features and configu- 
rations may occur as the dominant form of a phase transition. 

In small systems the system usually expands into the vacuum, and freezes out, thus the 
final state is out of thermal and mechanical equilibrium. Fluid dynamics cannot be applied 
at and after freeze-out and even in a short period before freeze-out, when the assumption of 
local thermal and mechanical equilibrium are not fully satisfied. Connections of freeze-out 
and hadronization are discussed in section |3|. 

In small finite systems among the configurations of possible instabilities we cannot 
neglect the ones associated by the outer surface of the system, which is of negligible im- 
portance for large macroscopic systems but may be dominant for small systems. 

2. Macroscopic Phase Transition Dynamics 

2.1. Slow dynamics - rapid phase transition 

When we calculate the speed of the phase transition proper, we have varying constraints 
and conditions. The speed of phase transition comes into question only if the dynamics of 
the evolution otherwise is so fast that it competes or exceeds the phase transition speed. 

If the external speed is slow, we have sufficient time to have a quasi-static process and 
reestablish phase equilibrium at every stage of the dynamics. This also means that all other 
equilibration processes are also completed as these require less time and less interaction 
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than phase transition dynamics. 

Thus, in the case of a "slow" external dynamics and rapid phase equilibration the 
matter is in complete equilibrium, including phase equilibrium, and the EoS of the mat- 
ter having a first order phase transition is given by the Maxwell construction: we have a 
fully developed mixed phase, and the phase abundances are given by the fluid dynamical 
evolution. No extra information on dynamical processes is needed. 

2.2. Deviation from phase equilibrium 

Even in moderately fast dynamical situations we have small deviations from the ideal and 
complete phase equilibrium (Maxwell construction). This deviation leads to some delay in 
the creation of the new phase leading to supercooling or superheating, and extra entropy 
production. 

For heavy ion reactions the first attempt to explicitly evaluate the phase transition speed 
of the homogeneous nucleation process is described in refs. [ [j], g]. The homogeneous 
nucleation mechanism describes correctly the initial phases of the phase transition, where 
the abundance of the newly created phase is still small, and when the phase transition 
process is the slowest. 

Here a couple of remarks are necessary. To form bubbles or phases of the new phase 
of supercritical size one needs to establish several requirements. Pressure and temperature 
balance should be reestablished among the phases and this requires to establish the phase 
boundary and transfer the needed energy and momentum to the new phase. We cannot 
relax the requirement of pressure and temperature equilibrium if both before and after the 
formation of the new phase we assume local equilibrium and so fluid dynamical evolution. 

2.3. Initial state 

As mentioned above the phase transition speed does not come into play for really large, 
slow stationary systems. If the dynamics is supersonic shock or detonation fronts will have 
the same initial and final parameters, and only the shock front profile may be affected by 
the phase transition dynamics. 

This is the typical type of approach also in heavy ion reactions up to a few hundred 
A ■ MeV colliding energy. 

3. Phase Transition Dynamics in Small Systems 

3.1. Final hadronization at freeze-out from QGP 

In principle it is possible that our system freezes out before kinetic and/or phase equilibrium 
is established. Then, we end up in a system of two or more phases, where none of them is 
equilibrated kinetically, so we have no partial pressures and temperatures. 

The case of freeze-out directly from QGP is very special, because hadronization must 
always be completed by the end of freeze-out, as QGP or free quarks and gluons never 
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reach the detectors. So, we must have completed phase transition to a single phase even if 
no thermal or mechanical equilibrium is established by the freeze-out. 

Recently final freeze-out and hadronization was discussed in a series of publications [ 
ILi' @> H improving essentially the previous Cooper-Frye freeze-out description 
[ [Tl] triggered by ref. [ |l2|| which suggested a solution to a long standing problem in 
freeze-out description, but did not offer a complete solution. 

In this series of works it is proven that in space-like (do^da^ = — 1) freeze-out the post 
freeze-out local momentum distribution cannot be a thermal distribution (because p^da^ > 
should be satisfied), and the earlier suggested cut-Juttner distribution is physically inad- 
equate. A simple kinetic model calculation provided an alternative non-equilibrium mo- 
mentum distribution, /(/?,e,n,w'',Ja'') (see Fig. |l]). Using this distribution it was shown 
that the freeze-out problem can be solved avoiding all previously mentioned problems. 




1 2 1 2 1 2 3 
p [100 MeV] 

Fig. 1. Post freeze-out momentum dis'tribution, ff ree (x : p), in the rest frame of the freeze- 
out front obtained in [ g]. Gradual kinetic freeze-out along x axis is described in a kinetic 
model. A, B and C correspond to freeze-out distances Q.2X, 3A-, 100A. respectively and 
u> rfg\ x = = (1 1 0-5, 0,0), where RFG - rest frame of gas, A, is scaling parameter of the order 
of mean free path. The numbers in the contours are in arbitrary units. The distribution 
is asymmetric and elongated in the freeze-out direction, x. This may lead to a large-/?, 
enhancement, compared to the usual Jiittner assumption used in most previous calculations 
as a freeze-out distribution. 

Most importantly the conservation laws must be exactly satisfied: 

[T fN da fl ]=0, [Woj=0, [S'Wa,,] > 0, (1) 

where the square bracket stands for the difference of the post and pre freeze-out values. 
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Neither these conditions nor the adequate post freeze-out non-equilibrium distributions 
were used in the earlier Cooper-Frye freeze-out models. If these conditions are satisfied 
at the freeze-out surface and the post freeze-out distribution, /(p,e,n,M'',^/o , ") is properly 
chosen the so called Cooper-Frye /ormM/a can still be used. 

The post freeze-out non-equilibrium distribution should nevertheless be evaluated in 
an adequate nonequilibrium dynamical model. The simple kinetic models used in refs. [ |J 
0- B 3 @ can ^ e considered as a first step only, which does not address important further 
details as post freeze-out flavor abundances, rapid and sequential freeze-out mechanisms, 
etc. A realistic and detailed freeze-out model requires a consistent dynamical model which 
generates such a distribution. 

Present string and parton cascade models if they can handle both the pre and post 
freeze-out phases realistically (!) may be used for this purpose with proper care. First 
attempts of coupling hadronic string models to the end of FD are performed already. [ |l3[ 
However, these models do not handle the phase transition dynamics, so satisfying the last 
condition, the requirement of entropy increase, is nontrivial in this model. 

Note that sudden or rapid hadronization from QGP with entropy increase usually pos- 
sible only in most model calculations if the plasma is sufficiently supercooled. 



3.2. Phase transition mechanisms 

The mechanisms of phase transitions have a large variety and consequently their dynamical 
features are also different. Experiments on strangeness and on two particle correlations 
suggest for some years by now, that hadronization and freeze-out is a rapid process, as 
the final observed system size is small and strange particle abundances are large. This 
contradicts to homogeneous processes in thermal and approximate phase equilibrium which 
were studied in the beginning of the 90s. 

In ref. [ ^| it was demonstrated that rapid QGP hadronization with small volume 
increase is possible even if we require energy- and momentum conservation and non- 
decreasing entropy. However, in ref. [ [|] no explanation was presented which would enable 
such a transition, which poses a problem as the earlier proposed homogeneous nucleation is 
not fast enough to support such a rapid transition. Alternative homogeneous processes (e.g. 
spinodal decomposition) were not studied in detail but as far as near thermal equilibrium is 
maintained it is difficult to imagine that these can support a qualitatively faster process. 

The need to find other processes which can support rapid hadronization was obvious 
by the mid 90s. The first attempt was to fully relax the requirement of thermalization, even 
the existence of temperature, and find a non-thermal, field theoretical, mechanism for the 
hadronization [ |]]. This connection gave then renewed activity in the study of fluctuations 
in field theories and in study of DCC. However, most of these studies were still considering 
a homogeneous transition. 



4. Phase Transition in the Initial State 



At highly ultra-relativistic energies, when QGP is formed, the pre collision initial state is 
far out of any thermal or mechanical equilibrium, and equilibrium can only be established 
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in the QGP phase where the number of degrees of freedom are sufficiently large and the 
interaction frequency is also large so that equilibration is expected to be established in a 
tenth of a fm / c. 

Nevertheless, the preceding dynamics cannot be treated by any thermal or fluid dynam- 
ical model, and we need QCD based effective models which originate from the observations 
gained from particle and heavy ion physics experiments at these energies. 

Just like connecting the hydrodynamic stage and freeze-out to each other on a 3 di- 
mensional hypersurface a detailed description of an energetic heavy ion reaction requires 
a Multi Module Model, where the different stages of the reaction are each described with 
suitable theoretical approaches. It is important that these Modules are coupled to each other 
correctly: on the interface, which is a 3 dimensional hyper-surface in space-time with nor- 
mal da^, all conservation laws should be satisfied and entropy should not decrease, eq. ([]]). 
These matching conditions were worked out and studied for the matching at FO in detail in 
refs. [g 10110. 

The initial stages are the most problematic. Frequently two or three fluid models are 
used to remedy the difficulties, and to model the process of QGP formation and thermaliza- 
tion. [ H Hi Here, the problem is transferred to the determination of drag-, friction- 
and transfer- terms among the fluid components, and a new problem is introduced with 
the (unjustified) use of EoS in each component in a nonequilibrated situations, where EoS 
does not exist. Strictly speaking this approach can only be justified for mixtures of nonin- 
teracting ideal gas components. Similarly, the use of transport theoretical approaches as- 
suming dilute gases with binary interactions is questionable, as due to the extreme Lorentz 
contraction in the CM. frame enormous particle and energy densities with the immediate 
formation of perturbative vacuum should be handled. Even in most parton cascade models 
these initial stages of the dynamics are just assumed in form of some initial condition, with 
little justification behind. 

All string models had to introduce new, energetic objects: string ropes [[l7[ |I8( ], quark 
clusters [ |l^j, fused strings [ |2(J], in order to describe the abundant formation of massive 
particles like strange antibaryons. Based on this, we describe the initial moments of the 
reaction in the framework of classical (or coherent) Yang-Mills theory, following ref. [ [H]] 
assuming larger field strength (string tension) than in ordinary hadron-hadron collisions. In 
addition we now satisfy all conservation laws exactly, while in ref. [ |TJ infinite projectile 
energy was assumed, and so, overall energy and momentum conservation was irrelevant. 

4.1. Coherent Yang-Mills model 

Our basic idea is to generalize the model developed in [ ^TJ, for collisions of two heavy 
ions and improve it by strictly satisfying conservation laws. First of all, we would create 
a grid in [x,y] plane (z - is the beam axes, [z,x] - is reaction plane). We will describe the 
nucleus-nucleus collision in terms of steak-by-streak collisions, corresponding to the same 
transverse coordinates, {x,-, yj}. We assume that baryon recoil for both target and projectile 
arise from the acceleration of partons in an effective field F MV , produced in the interaction. 
Of course, the physical picture behind this model should be based on chromoelectric flux 
tube or string models, but for our purpose we consider F^ v as an effective abelian field. 
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Phenomenological parameters describing this field must be fixed from comparison with 
experimental data. 

Let describe the streak-streak collision. 

( 2 ) 

i i 

3„£X = 0, i = l,2, (3) 

is the baryon current of /th nucleus (we are working in the Center of Rapidity Frame 
(CRF), which is the same for all streaks. The concept of using target and projectile refer- 
ence frames has no advantage any more). We will use the parameterization: 

nf = Piul , ul = (coshy,-, sinhy;) . (4) 

is a energy-momentum flux tensor. It consists of five parts, corresponding to both 
nuclei and free field energy Tp[ (also divided into two parts) and one defining the QGP 
perturbative vacuum. 



Bg» v , i=l,2, (5) 



where B - is the bag constant, the equation of state is P,- = c^e,, where e, and P, are energy 
density and pressure of QGP. In complete analogy to electro-magnetic field 

F? v = 3 V A? - VA] = ( ° ~ ai ) , where o, = 3 3 A° - 3°A? , (6) 



a, 

7>, ( >v = -guvXF,f + y* f F ' l a s 9v^f , where = —F^Ff . (7) 



The string tensions, a,, have the same absolute value a and opposite sign (in complete 
analogy to the usual string with two ends moving in opposite directions), and o, will be 
constant in the space-time region after string creation and before string decay. 

We received analytic solutions of the above equations using light cone variables [53], 
x ± = t ± z. Following [ |23|] we assume that target variables, e\, y\ , pi , are functions 
of x~ only and projectile variables of x + only. At the time of first touch of two streaks, 
t = 0, there is no string tension built up yet. We assume that strings are created, i.e. the 
sting tension achieves the value a at time t = to, corresponding to complete penetration of 
streaks through each other. 

4.1.1. Conservation laws — string rope creation 

In light cone variables the baryon current conservation, eq. (|j) may be rewritten as 

d-n^+d+n^ = 0. (8) 
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So, we have a sum of two terms, depending on different independent variables, and the 
solution can be found in the following way: 3_n^ = a, 8 + nJ = —a , with = ax~ + 
(«i)o, « 2 — —ax + + («2)o • Since both and are positive (and also more or less 
symmetric) we can conclude that for our case a = 0. Finally: =p\e~~ yx — po<?- v ° , n 2 = 
p 2 e>'2 = p oe >'o j an d pi = p oe yo+vi ( p 2 = po^'o-w . 

As mentioned before, after string creation, i.e. t > to, and before string decay we 
choose the string tensions in the form: a 2 = — CTi = O > . With this choice and the 
Lorentz gauge condition we take the vector potentials in the following form [ [22| ^3|] : 

A+ =0, A~ = -2ax" ; A+ = -2ox + , A 2 = , (9) 

where we used the parameterization: 

O=A(^) 2 p y/hk, (10) 

\ m / 

where l\ and h are the initial streak lengths. We are working in the system, where h — c= 1 . 
The typical values of dimensionless parameter A are around 0.045 — 0.055. The above 
parameterization is arbitrary in the sense that the requirements of the right dimension and 
grid size independence do not completely fix it. Above parameterization has been checked 
to work in the energy range £0 = 10 — 3000 GeV per nucleon. Notice, that there is only 
one free parameter in parameterization (|Io|). The typical values of a are 5 — 12 GeV / fm 
for Eo = 100 GeV per nucleon. These values are consistent with the energy density in 
non-hadronized strings, or "latent energy density" which is on the average 8 GeV/fm 3 .[ 

As eq. (g) has a source term we do not know what the really conserved quantities are. 
Using the solution for n[ and , we can define new energy-momentum tensor T^ v , such 
that 

d„J* v = 0, T " v = £ ?T + T%x = I (7f - AX + g" V Afn ia ) +Bg>» (11) 

i i 

Now the new conserved quantities are 
Qo = j f m dV = AxAy^Jj^dz , Q3 = J f m dV - AxAyJ^Jj^dz , (12) 

where Ax Ay is the cross section of the streaks. 

Then the trajectories of nucleons (or cell elements) for both nuclei, the energy density 
and baryon density distributions can be obtained analytically and are given in [ ^3|. 

4.1.2. Recreation of the matter 

As we may see from the trajectories nucleons (or cell domains) will keep going in the 
initial direction up to the time t = tjj ur „, then they will turn and go backwards until the two 
streaks again penetrate through each other and new oscillation will start. Such a motion is 
analogous to the "Yo-Yo" motion in the string models. Of course, it is difficult to believe 
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that such a process would really happen in heavy ion collisions, because of string decays, 
string-string interactions, interaction between streaks and other reasons, which are quite 
difficult to take into account. To be realistic we should stop the motion at some moment 
before the projectile and target cross again. 

We assume that the final result of collisions of two streaks after stopping the string's ex- 
pansion and after its decay, is one streak of the length If with homogeneous energy density 
distribution, et, and baryon charge distribution, Pf, moving like one object with rapidity 
yf. We assume that this is due to string-string interactions and string decays. As it was 
mentioned above the typical values of the string tension, o, are of the order of 10 GeV / fm, 
and these may be treated as several parallel strings. The string-string interaction will pro- 
duce a kind of "string rope" between our two streaks, which is responsible for final energy 
density and baryon charge homogeneous distributions. Now it is worth to mention that 
decay of our "string rope" does not allow charges to remain at the ends of the final streak, 
as it would be if we assume full transparency. 

The homogeneous distributions are the simplest assumptions, which may be modified 
based on experimental data. Its advantage is a simple expression for e/, Pf, y/. 

The final energy density, baryon density and rapidity, ef,pf and yf, should be deter- 
mined from conservation laws. Unfortunately, the assumptions we made above oversim- 
plify the real situations and do not allow us to satisfy exactly all the conservation laws. The 
reason for this is well known and has been discussed in the Refs. [ |[ 0, f| |[ |i0[ ]: two 
possible definitions of the flow, Eckart's and Landau's definitions, are not always identical. 
The exact conservation of the energy and momentum gives for the final rapidity: 

(M 2 (l +4) +2cy Q ) + J(Ml(l +c 2) +2 C 2 M 2 ) 2 +4c 4 M 2 (M 2 _„2 } 

cosh yf - L = 2 { i + 4W-4) ' (13) 

where we neglected B A If next to Qq/ Ax Ay and introduced the notation M = (lj + 
h)/(h — h), uq = tanhyo is the initial velocity. The exact conservation of the baryon four- 
current gives: 

uq j M 2 

tanh y/ , £ =-,^ C osh y/!£ = ^— j. (14) 

It is interesting that if we put Cq = the eq. ( |j"3| ) becomes identical to eq. (0). For more 
details see [ ^3J . 

We follow the Refs. [ ^2|, ^3[], where the yj = yfx has been chosen. In this case the 
expressions for the e f and p t are: 

Qo 



e f = ^ (15) 

f ((l+cl)coshiy f -cl)l f > 

Pojh+h) 

Pf = -, r — ■ (lo) 

Ifcoshyf 

So, the streaks move until they reach the rapidity y, = yf. Later the final streak starts 
to move like one object with rapidity yf. 
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4.2. Initial conditions for hydrodynamical calculations 

We are interested in the shape of QGP formed, when string expansions stop and their matter 
is locally equilibrated. This will be the initial state for further hydrodynamical calculations. 
We may see in Figs. ^, that QGP forms a tilted disk for b ^ 0. So, the direction of fastest 
expansion, the same as largest pressure gradient, will be in the reaction plane, but will 
deviate from both the beam axis and the usual transverse flow direction. So, the new flow 
component, called "antiflow" or "third flow component", may appear in addition to the 
usual transverse flow component in the reaction plane. With increasing beam energy the 
usual transverse flow is getting weaker, while this new flow component is strengthened. 
The mutual effect of the usual directed transverse flow and this new "antiflow" or "third 
flow component" leads to an enhanced emission in the reaction plane. This was actually 
observed and studies earlier. One should also mention that both the standard transverse 
flow and new "antiflow" contribute to the "elliptic flow". 
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Fig. 2. The Au+Au collisions, e = 100 GeV/nucl, b = Q.5{RV^R2j,'A l = 0.05 (parameter 
A introduced in (10)), y = (ZX plane through the centers of nuclei). We would like to 
notice that final shape of QGP volume is a tilted disk w 45°, and the direction of the fastest 
expansion will deviate from both the beam axis and the usual transverse flow direction, and 
might be a reason for the third flow component, as argued in [ ETI] . 



5. Critical Fluctuations 

Fluid Dynamics inherently describes the average behavior out the thermodynamical en- 
semble characterizing a system. Event by event fluctuations can nevertheless be included 
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in FD, and one can generate an ensemble of FD events. For mesoscopic systems the fluctu- 
ations are not negligible, and near the critical point these can even dominate the dynamics. [ 
p8[ |2^, [30|] Critical fluctuations may signal the vicinity of the critical point in a phase 
transition, thus may serve as a QGP signal. 

One may use the Dissipation-Fluctuation Theorem [ ^9|] to generate fluctuating Langevin 
forces on physical ground. Note that although this problem for the case of Navier-Stokes 
equation was addressed already by Landau in 1957, a practically usable approach was only 
worked out recently. 

In heavy ion reactions the latent heat of the phase transition released at the final 
hadronization may feed fluctuations, among other possibilities also critical fluctuations in 
FD. The first methodological steps were made in this direction, but a realistic and experi- 
mentally verifiable prediction of observable fluctuations is still some way ahead. 



6. Conclusions 



We recalled arguments for a rapid freeze-out in heavy ion collisions, which coincides with 
the hadronization of QGP, and presented a considerably improved method for the idealized 
description of the freeze-out. 

Based on earlier Coherent Yang-Mills field theoretical models, and introducing effec- 
tive string tension parameters based on Monte-Carlo string cascade and parton cascade 
model results, a simple model is introduced to describe the pre fluid dynamical stages of 
heavy ion collisions at the highest SPS energies and above. The model predicts limited 
transparency for massive heavy ions, as a consequence of collective effects related to QGP 
formation. These collective effects in central and semi central collisions lead to much less 
transparency than earlier estimates. The resulting initial locally equilibrated state of matter 
in semi central collisions takes a rather unusual form, which can be then identified by the 
asymmetry of the caused collective flow. Our prediction is that this special initial state may 
be the cause of the recently predicted "antiflow" or "third flow component". 
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